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We study the contribution of spectator massive scalar fields to the inflationary density pertur- 
bations through the universal gravitational coupling. We find that such contribution has several 
remarkable properties: it does not decrease as the mass of the spectator field increases; it has a 
significant size and cannot be turned off by any adjustable parameters; and it applies to all mas- 
sive scalars existed during infiation, making the overall effect unexpectedly large. As a result, the 
primordial density perturbations are anomalously sensitive to the high energy physics. 



Fields with masses much larger than the Hubble pa- 
rameter are abundant in any realistic models of inflation. 
The usual description of inflation models are given in 
terms of one or several light degrees of freedom in terms 
of the low energy effective field theory, such as done since 
the earliest models [HE]- The implicit assumption is that 
the heavier fields are effectively decoupled, analogous to 
the situation in particle physics. 

For example, although the massive fields can be cou- 
pled to the inflaton through some model-dependent direct 
couplings in the matter Lagrangian, it has been shown 
through examples that their contribution to inflationary 
density perturbations is inversely proportional to some 
power of mass [3l E] . Therefore heavier flelds are less im- 
portant, consistent with the intuition. Of course, certain 
low-lying massive flelds may become more accessible if 
the model contains sharp features [H [3] or special ini- 
tial conditions [7] that introduce higher energy scales, 
analogous to the injection of energies in particle physics 
colliders. But this does not affect the conclusion for the 
much heavier flelds. 

While the direct couplings may be present rather 
generically and exhibit many interesting phenomenolo- 
gies, they are not the minimal case. Namely, their 
strength is parameterized by certain coupling parame- 
ters, and there always exists a model-building limit in 
which we can turn them off. However, for inflaton after 
we turned off all the direct couplings, there is still one 
coupling remaining - the gravitational coupling - which 
couples all fields existed during inflation to the inflaton 
and forms the baseline for all inflationary models. This is 
the subject of this paper. We would like to compute the 
contribution of the massive flelds to the inflationary den- 
sity perturbations through the gravitational coupling. At 
one-loop level, this is given by the two diagrams in Fig. [T] 
Due to the special properties of gravity, we will see many 
dramatic differences from the other cases. 

We start with one massive scalar a with mass M in a 
slow-roll inflation model. The matter Lagrangian is 



(a) 



(b) 



FIG. 1: Two one-loop diagrams. C, is the scalar perturbation, 
a is the massive field. 



The (j) is the inflaton and Vsy{4>) is the slow- roll potential, 
with slow-roll parameter e = —H/H^. The massive fleld 
a is always decoupled from the inflaton in the matter 
Lagrangian and classically stays at ao(t) = 0. 

To study the density perturbations in the in-in formal- 
ism [3 IH], we use the ADM formalism [TIIIITT] . 

ds^ = -N^dP 4- h,j{dx' -f N'dt){dx^ + Wdi) . (2) 

We would like to write the perturbation theory in terms 
of the uniform inflaton gauge (C-gauge) 

% = a2(f)e2«(*-*)j,, , 0(t,x) = 4>o{i) , a - a(£,x) , (3) 

because in this gauge the physical variables are most 
transparent. For single field inflation, C, is the scalar per- 
turbation that conserves at the supcrhorizon scales and 
directly corresponds to observables in the Cosmic Mi- 
crowave background and Large Scale Structures, fn our 
case the massive field a \s a spectator and does not play 
any role in reheating. Therefore the physical meaning of 
C remains the same after reheating, although as we will 
see it is no longer conserved outside the horizon during 
inflation. The a in this gauge hence becomes the physical 
massive fleld that we will integrate over in the loops. 

The simpler way to see the relevant terms in the cu- 
bic and quartic perturbative Lagrangians, especially the 
orders of magnitude of these terms, is however not to di- 
rectly work in this gauge. In the ^-gauge, at the quartic 
order there are a large number of terms whose mutual 
cancelations are obscured by many integrations by part. 
Instead, it is easier to first work in the spatially flat gauge 
((50-gauge), 

= a^{t)5ij , (j) = Mt) + H{t,^) , cr = a{t,x) , (4) 

and then transform it to the ^-gauge through the coor- 
dinate transformation. 



(1) 



i = t + St{t,x') , x' = x' + 5x\t,x') . 



(5) 



See [12] for details. We have distinguished the two dif- 
ferent gauges by adding hats to the variables and coor- 
dinates in the i^-gauge. The transformation between the 
fields in the two gauges are given by [TH] 

6<j){t, x') = Mt){St(i) + 5t(2) + 5^^)) 

+ T^M^Ai) + 2'5i(i)<5t(2)) + J , (6) 
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(7) 

where the subscripts in St and Sx' label the perturbative 
orders, and we have listed terms up to the third order. 
For the purpose of this paper, the most important term 
is 

^t(i) = -C/H . (8) 
The orders of magnitude of the other terms are 

St^2) = O(C'), ^i(3) = 0{C\ ^'0, Sx(2) = OiC). (9) 

In the slow-roll expansion, the leading order La- 
grangian in the (5(/)-gauge is particularly simple, 

(10) 

The cubic and quartic Lagrangians are of order 
0{y^6(j)a'^) and 0{S(p'^a'^), respectively. So they are e- 
suppressed in terms of the C-gauge. 

We now apply the gauge transformation. Higher or- 
der terms generated by ^ are suppressed by e. So the 
leading terms in £3 and £4 in the ^-gauge are entirely 
generated by Among all these terms, terms of order 
0{(a^) generated from £2 vanish except for temporal 
total derivative terms, because they are proportional to 
the equation of motion after integration by part. Same 
for terms of order 0{('^a'^) generated by the third order 
terms in ^ . But the quadratic order terms in ^ , acting 
on all cr's in ( [To| ), can generate terms of order C(C^(T^). 
To summarize, up to temporal total derivative terms |12j . 
the leading terms in the C-gauge Lagrangian are 



£3 = OieCa') , 
„3 r 



(11) 

(12) 
(13) 



Here and below we drop all the hats for simplicity. Each 
derivative on a produces a factor of momentum cutoff 
A in loops. So £4 is apparently of order A"*. However, 



using spatial integrations by part and the temporal to- 
tal derivative terms, we can reduce it to order A^. This 
avoids inconvenient cancelations between large numbers 
in the calculation. Finally we perform a Legendre trans- 
form and obtain the following kinematic and leading or- 
der interaction Hamiltonian, 



Ho - ea'c' + eaidCf + + ^(a.a)^ + "^M^a^ 



(14) 
(15) 



We have only listed three terms in Hi, which will be used 
in this paper as the representative examples. A full list 
of terms is presented in [12,. The dots in (15) include 



many other similar quartic order terms from (131 and 



those generated by the Legendre transform due to £3. 
The dots also include the leading cubic terms related to 
the temporal total derivative terms in £3, and the precise 
treatment of such terms can be found in 1121 . 



We now use the three representative terms in ( 15 ) 



to compute the contribution of the massive field to the 
power spectrum. 

The quantization of the fields and the mode functions 



followed from the kinematic Hamiltonian ( 14 1 are the 



same as in Ref. |4j. In terms of the Fourier components, 
we quantize 

Ck Mkflk + ulk^-k ' = Wk^k + ^'-k^-k ' (16) 

where ak and 6k are independent of each other and each 
satisfy the usual commutation relations. The mode func- 
tions satisfy the equations of motion 



-Mk -f /c^Mk — , 



(17) 
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v'k - + k^vi, + -jj^Vk , (18) 

where r is the conformal time, dt = adr, and the primes 
denote derivatives with respect to r. The solutions are 



H 



Mk 



:(1 -|- ikr) e 



-ikr 



(19) 



where ^ = -^/M^T^^^'^-W^- Both of the mode functions 
approach the Bunch-Davies vacuum in the short wave- 
length limit. The power spectrum is defined as 



(CkCk') = §(2^f^'(k + k') 



At the tree level, 



(21) 



(22) 
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At the CMB scales ko = 0.002Mpc-\ w 6.7 x 10"^ 

Using the in-in formalism, the diagram (b) in Fig. [l] 
gives the following correction to the tree-level two-point 
correlation function (C^), 

A{C) D 2i 04(Te„d))' {27r)'S^k + k') 



dr 



oo 
„2 



3a 



H 
(27r)3 



(27r)3' 



The integration over the loop momentum p is power- 
law divergent at UV {p — > oo) and needs to be renormal- 
ized. As in particle physics, we impose a UV cutoff scale 
in terms of the physical momentum A, which in terms 
of the comoving momentum is pA = — A/(ifr). This is 
the energy scale below which our effective field theory 
description is valid for both the inflaton and the massive 
field. The detailed value of A depends on how the field 
theory is UV completed and will be model-dependent, 
and we use A to parameterize our ignorance. 

The integration over the conformal time r is logarith- 
mically divergent at IR (rond 0). This divergence has 
the physical interpretation that the scalar perturbation 
is no longer constant after the Hubble horizon exit. It re- 
ceives constant contributions from the massive field loop 
per Hubble time. There exists a natural IR cutoff to 
this divergence - the end of inflation. The result there- 
fore is proportional to ln(— 2fcTend)) which is the number 
of efolds Ne till the end of inflation after the mode 2k 
crosses the horizon. 



The integration (23) can be done exactly and the re- 
sults are given by special functions |12j . The most im- 
portant features and physical meaning of the results, 
however, are best illustrated by looking at some sim- 
plification limits. In the large momentum cutoff limit, 
A/H ^ (M/H)'^, we can approximate the Hankel func- 
tion in Vp by its UV expansion. 



g-j7r/4g7r/j/2 



irpr 



(24) 



From ( 23 1 we get 



A4 



Pc 



19 3M2 
Sn^H^M^^e \ 12 ^ 



ln(-2fcre„d) . (25) 



It is also interesting to look at the limit A/H ^ M/ H 
in which we can approximate the Hankel function by its 
IR expansion. 



where 5 is an unimportant phase. In this limit. 



(26) 



A4 



( 



3M 



A 



Sn^H^M^^e \ A 15M 



ln(-2A;Te„d) ■ (27) 



a: 
<1 
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(23) 

FIG. 2: Correction (|23l) to the curvature perturbation power 



spectrum as a function of isocurvaton mass, with A — 10 H. 



For A ^ M, the actual behavior is similar in order of 
magnitude to both expressions (25) and (27). In Fig. [2] 
we compare the exact integration with these two approx- 
imations. The other terms not listed in ( |15[ ) give similar 
behaviors with different numerical factorsfl^ . 

These results have several remarkable properties. 

Firstly, comparing to the direct coupling case [21 H] in 
which the corrections are inversely proportional to power 
of A/, here I^Pq/Pq is either approximately independent 
of M or increases as M increases. This is because, in 
the interaction Hamiltonian, the gravitational coupling is 
proportional to the mass-squared or momentum-squared 
of the massive field. At the same time, the oscillatory 
factor of the massive field cancel between the particle 
and anti-particle in the loop. Such a behavior is in sharp 
contradiction to the usual intuition that the massive field 
and high energy physics are decoupled from the low en- 
ergy effective theory. 

Secondly, again different from the direct coupling case, 
here there is no adjustable parameters one can use to turn 
off the coupling. Despite the Planck-mass suppression, 
the correction to the power spectrum is actually quite 
significant even if we consider single massive field. To 



see this, we can write (25) as 



AP, 



C 



(28) 



where we have used ( p2| ) and ln(— 2A:Tend) ~ N^. ^ 60. 
For as low as A > M > O(10^)i/, the loop contribution 
from the massive field already starts to exceed the tree- 
level result, and the usual perturbation theory breaks 
down. An even stronger constraint could be obtained by 
examining the contribution of the logarithmic term to 
the tilt of the power spectrum. 

Lastly, and more importantly, the results apply to all 
scalar fields existed during inflation [14 . Among them 
the most interesting ones for our purpose are the massive 
flelds with M > H. The existence of such fields are 
guaranteed in any realistic models. Furthermore, they 
often show up as towers of states in string theory models, 
with the number of states growing with the masses. 
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To summarize, the inflationary density perturbations 
turn out to be anomalously sensitive to the high energy 
physics, through the gravitational coupling between the 
inflaton and all the other states. 

To further illustrate these unusual features, we apply 
the results to an example of extra-dimensional models. 
We consider D-dimensional spacetime compactified on a 
{D — 4)-dimensional compact space, such as a (£> — 4)- 
torus. The sizes of the compact dimensions are of order 
L, and the fundamental Planck mass is M^. The 4D 
Planck mass is given by 



(29) 



The extra dimensions give rise to towers of massive states 
such as the scalar Kaluza-Klein states. The mass is given 
by the KK number n^, 



1 

Z2 



_D-4 



(30) 



The KK states can be described by the 4D effective 
field theory until the mass of the highest KK states 
becomes of order Mjj, namely the KK numbers reach 
\/ nf — rimax ~ MuL. Any inflationary models in this 
setup introduce the gravitational coupling between these 
KK states and the inflaton. The momentum cutoff in 
the loop is A ^ Mn- Because the density of the KK 
states grows with the mass, the dominant contributions 
come from the highest states with mass of order the cutoff 
scale. The total number of such states is 



D-4 



[MdL) 



Ml 



The net contribution of all these states is 



C 



MlN, 



(31) 



(32) 



Because Md > H, it is clear that this quantity is al- 
ways much larger than one. Therefore the conventional 
tree-level result is by no means the obvious leading order 
result, and it is essential to reconsider the calculation 
by taking into account more details of the high energy 
physics. 

There are many important issues remaining to be in- 
vestigated |15| . For example, it is important to under- 
stand how the density perturbation computation should 
work in cases where the first order gravity-mediated loop 
contribution becomes more important than the tree level 
result. One possibility is that one may be able to resum 
all the loop diagrams. This will completely change the 
prediction in the non-perturbation regime. Another pos- 
sibility is to introduce supersymmetry so that the loop 
contribution from fermionic massive states exactly can- 
cel the bosonic ones above the supersymmetry breaking 



scale, similar to what may happen to the Higgs mass hi- 
erarchy, or the cosmological constant problem. It would 
be interesting if supersymmetry were essential for mak- 
ing sense of such basic computations of density pertur- 
bations. It is important to extend the scalar in the loop 
to other type of fields. It is also important to consider 
how the calculations of the non-Gaussianities and tensor 
power spectrum, and how other inflation models besides 
the slow-roll models, are affected by the effects we con- 
sidered here. 
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